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Abstract. We investigate one pion production processes within the Giessen Boltzmann-Uehling- 
Uhlenbeck (GiBUU) coupled channel transport model. Our calculations for integrated and differ- 
ential cross sections for realistic experimental neutrino fluxes are compared to the data recently 
Cn ■ provided by the MiniBooNE collaboration. 
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INTRODUCTION 



Recently the MiniBooNE and K2K collaborations have published data on charged [1] 
and neutral [2] pion production in CC neutrino scattering. For the first time the pion 
energy and angle observables were presented. Generally, these experiments report cross 
sections which are noticeably higher than those expected from any conventional theo- 
retical approach [3]. 

In this paper we investigate by detailed comparisons with data on one pion production 
how far a realistic impulse approximation model of these processes can go in explaining 
all the available data. For this purpose we employ the GiBUU model developed as a 
transport model for nucleon-, nucleus-, pion-, and electron- and neutrino-induced colli- 
^ ! sions from hundreds MeV up to tens GeV. Thus, the model allows one to study various 

types of processes on nuclei within a unified description. This is particularly important 
for broad-band neutrino experiments which inherently sum over many different reaction 
types. The code is written in modular FORTRAN and is available for download as open 
source [4]. Our results for QE scattering and NC pion production are presented in [5]. 



GIBUU TRANSPORT MODEL 

The MiniBooNE energy flux [6] has an average neutrino energy 0.7 GeV, with the flux 
becoming very small above 2 GeV. Several different channels are important for de- 
scribing the neutrino interactions at these energies: quasi-elastic scattering, resonance 
production (with the A [Pt,j,(1232)] giving the biggest contribution) and non-resonant 
one pion production. All of them are described within the GiBUU model, the de- 
tails are given in [7]. Our elementary input is described in [8]. We emphasize, that 
in describing QE scattering the world average value for the nucleon the axial mass 
Ma = 0.999 GeV [9] is used; in describing resonance production and background pro- 
cesses parameters were tuned to the ANL data. 



Nuclear targets are modeled as follows. The struck nucleus is considered as a collec- 
tion of off-shell nucleons. Each nucleon is bound in a mean-field potential, which on 
average describes the many-body interactions with the other nucleons. The phase space 
density of nucleons is treated within a local Thomas-Fermi approximation. In the im- 
pulse approximation, an incoming lepton interacts with a single bound nucleon, with the 
interaction vertex being the same as in the case of a free nucleon. Thus, the calculations 
do not contain any so-called 2p-2h interactions [10]. 

In the GiBUU code it is also possible to take into account the nuclear medium 
corrections to the resonance widths and spectral functions. In this paper, however, they 
are neglected. Final State Interactions (FSI) are implemented by solving the semi- 
classical Boltzmann-Uehling-Uhlenbeck (BUU) equation. It describes the dynamical 
evolution of the phase space density for each particle species under the influence of the 
mean field potential, introduced in the description of the initial nucleus state. Equations 
for various particle species are coupled through this mean field and also through the 
collision term. This term explicitly accounts for changes in the phase space density 
caused by elastic and inelastic collisions between particles. 

FSI decrease the cross sections as well as significantly modify the shapes of the 
final particle spectra. Such modification was experimentally observed in photo-pion 
production [11, 12, 13] and is described by the GiBUU with a high accuracy. 

For a detailed review of the GiBUU model see [7]. 

ONE PION PRODUCTION AND THE ORIGIN OF PIONS. 

Here we present our calculations on a CH2 target and compare the results with the data 
from the MiniBooNE experiment [1, 2], Integrated cross sections versus neutrino en- 
ergy for the charged current ln + and in production are shown in Fig. 1. As in the 
MiniBooNE experiment, the one pion events are defined as "observable one pion pro- 
duction", i.e. events with one pion of a given charge and no other pions in the final state, 
regardless of which particles were produced in the initial neutrino vertex. 

The left upper panel in Fig. 1 shows the results for n + production. Comparison 
of the curves with and without FSI shows that the FSI do not change the energy- 
dependence of the curves. The curve with FSI (solid line) lies 20% below the curve 
without FSI (dash-dotted line). This reduction is mainly caused by the A absorption 
through NA — > NN scattering. Charge exchange process K + n — > n°p also depletes the 
dominant \k + channel. The side feeding from the reverse process gives minor relative 
contribution, because the initial n° production cross section is around 5 times lower. For 
the l7T° cross section in the left lower panel the same side feeding processes increase 
the cross section. This is, however, compensated by absorption, charge exchange to the 
K~ channel through n°n — >■ %~ p and other channels such as kN — > ZK,AK. The overall 
effect is nearly the same cross section with and without FSI. 

The right upper panel shows the origin of the ln + events. The most of them (dash- 
dotted line) come form the initial Delta resonance production and its following decay. 
Some events (dotted line) are background ones. A minor amount comes from the initial 
QE vertex (long-dashed line), which is only possible due to FSI, when the outgoing 
proton is rescattered. Here the main contribution is from the pN -> N'A — > N'N % 
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FIGURE 1. (Color online) Integrated cross section for MiniBooNE l7T + (upper panel) and l7T°(lower 
panel) CC production versus neutrino energy. Data are from [1, 2], The panels on the right show the 
composition of the calculated cross sections. 



reaction. The right lower panel shows the origin of the lft events. Here the background 
processes and the FSI play an even bigger role. 

We now turn to a discussion of a comparison with experiment. Already the curves 
without any FSI lie considerably (rj 25%) below the data, those with FSI included (solid 
line) are nearly by a factor of about 1 .5 (at 1 GeV) below the experimental data; at higher 
neutrino energies these discrepancies become even larger since the experimental cross 
sections rise more steeply with energy than the calculated pion cross sections. At 2 GeV 
the discrepancy between the results with FSI and the data amounts to a factor of ~ 2! For 
In production one observes the similar result: the data are considerably higher than our 
calculated values. The slope of the curve, however, is in agreement with the experiment, 
as one can conclude from comparison of the "with FSI" curve multiplied by a factor of 
1.5 (long-dashed curve) with the data. 

An observable that is approximately free of the uncertainties of the neutrino flux is the 
ratio of pion production to QE-like scattering at fixed neutrino energy. Here, of course, 
an additional uncertainty enters because the energies for QE and for one pion production 
are being reconstructed independently and from different underlying assumptions. The 
errors in the two independent energy-determinations thus add. 

This is essential since our calculations for the ratio of \n + to QE-like cross sections 
shown in the left panel of Fig. 2 gives a reasonably good agreement with the data up 
to about 1 GeV neutrino energy and even above that the discrepancy amounts to only 
about 10%. This uncertainty is well within the limits of the combined errors of flux 
and energy reconstruction. Note here that the experimental QE-like cross section most 
probably contains contributions from 2p-2h excitations while the theoretical does not. 
Including those in the calculations would lower the ratio by about 25% at 1 GeV. 




o 

O 0.2 



0.6 0.8 1 1.2 1.4 1.6 1.6 
E v , GeV 




0.6 0.8 1 1.2 1.4 1.6 1. 
E v , GeV 



FIGURE 2. (Color online) (Left panel) Observed single l7T + /QE-like cross sections ratio for CC 
neutrino scattering versus neutrino energy. Data are from [14]. (Right panel) Integrated cross section 
for MiniBooNE n~ CC production versus neutrino energy. 



Our prediction for the %~ cross section are shown in the right panel of Fig. 2. Without 
FSI the production of %~ is negligible, a few events may only come from higher mass 
resonances in the processes like R + — y pp° — >■ pll + il~ orR + — > pA° — > ppn~ (the curves 
are not shown because they would be indistinguishable from zero). Cross sections for 
both \%~ (solid curve), which is defined as one %~ and no other pions, and multi-7T~, 
which is defined as at least one n~ and any number of pions of other charges, are shown. 
The cross sections for %~ production appear to be around 20 times smaller than those 
for ln + . However, with the current statistics of the MiniBooNE experiment they should 
be measurable. 



LEPTON SPECTRA 



A more detailed information on pion production processes is contained in differential 
distributions. In Figs. 3, 4 we compare the recent MiniBooNE distributions versus muon 
kinetic energy, do/dT^, and squared four- momentum transfer, do/dQ 2 [1] with our 
calculations averaged over the MiniBooNE flux. 

As could be expected from the comparison with the integrated cross sections, our 
calculations with FSI (solid linse) and even without FSI (dash-dotted line) are also 
significantly lower than the experimental data for both ln + and \ir events. As in the 
case of the integrated cross sections, our calculations even without FSI, shown as dashed- 
dotted lines, are about 20% lower than the data. These results represent essentially only 
the free cross section folded with the Fermi momentum-distribution. The disagreement 
is therefore, not easy to understand. 

To make a shape only comparison for ln + , we multiply our "with FSI" calculations 
by a factor 1.6 (long-dashed line) and the "without FSI" by 1.2 (dotted line). These 
curves practically coincide, which means, that FSI hardly change the shapes of the 
distributions. The shape is in good agreement with the shape of the experimental data 
for the Q 2 distribution, and in reasonable agreement for the kinetic energy distribution. 
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FIGURE 3. (Color online) Distribution of the outgoing muons in their kinetic energy for MiniBooNE 
CC l7T + (left panel) and In (right panel) neutrino production. Data are from [1, 2]. 
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FIGURE 4. (Color online) Q 2 distribution for MiniBooNE CC 1 7T+(left panel) and l7T° (right panel) 
production. Data are from [1,2]. 



Recall here that for the integrated cross section we had a noticeable shape disagreement 
in this channel. 

For the l7T° channel we we multiply our "with FSI" calculations by a factor 2.0 (long- 
dashed line). Making a shape-only comparison, one observes noticeable deviations for 
both differential cross sections, despite a good shape-only agreement for the integrated 
cross sections. Our curve for the 7^ distribution is noticeably flatter than the data. The 
calculated Q 2 distribution is steeper. Since both low 7^ and low Q 2 originate from low- 
energy muons, the disagreement is hard to explain. 

CHANGE OF THE PION SPECTRA DUE TO FSI 



Fig. 5 shows our calculations for the kinetic energy distribution of the outgoing pions. 
As for the other distributions, our calculations with FSI are lower than the experimental 
data by a factor of 1.6 for n + and a factor of 2 for n . 

A distinctive feature of our calculations is the conclusion, that FSI should significantly 
change the pion spectra. The significant lowering of the cross section for T K + > 0. 12 GeV 
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FIGURE 5. Distribution of (left panel) the outgoing n + in their kinetic energy; (right panel) the 
outgoing 7T° in their absolute value of the 3-momentum for the MiniBooNE CC neutrino scattering. Data 
are from [1, 2]. 



and for p\ > 0.25 GeV, is a direct consequence of the pion absorption through tzN — » A 
with the following AN — > NN. Pion elastic scattering in the FSI also decreases the 
pion energy, thus depleting spectra at high energies and accumulating strength at low 
energies. For n° production an additional increase of the cross section at lower energies 
comes from the side feeding from the dominant n + channel as discussed above. The 
change of the shape of the spectra is similar to that calculated for neutral current In 
production in [15]. 

The predicted low-energy peak in pion spectra is also similar to that observed exper- 
imentally in (y, n ) production on nuclear targets [11, 13]. Its absence in the neutrino 
data is, therefore, hard to understand. At very low energies our curves rise as steeply as 
the data. After the peak is reached, however, our calculations predict a rather steep fall- 
off, while the data show only a moderate decrease. In this view, it would be extremely 
important to compare our calculations with the coming results from the CLAS experi- 
ment for pion distributions in electron scattering on nuclei (see the talk of S. Manly at 
these conference). 

ANY REASONABLE EXPLANATION FOR DISCREPANCIES? 



One possible explanation for the discrepancies observed could be too low elementary 
cross sections since the ANL and BNL data for pion production already differ, with the 
BNL data being higher by about 30% (recall that our calculations use form factors for 
the A^A transition that were fitted to the ANL data). However, even such an increase 
would not be sufficient to explain the large discrepancy by 60% and more. 

Another source of uncertainties is the determination of the neutrino energy which is 
done using quasi-free kinematics for on-shell Delta production on a nucleon at rest. This, 
however, have no influence of flux-averaged observables. One more point is that the ex- 
perimental data are presented after being corrected for the finite detection thresholds for 
muons, pions and nucleons. This may introduce additional dependence on a particular 
neutrino event generator used in a given experiment [16]. 



Another possible explanation for the discrepancy is an uncertainty about the actual re- 
action mechanism. Recently, several authors have raised the possibility - in connection 
with QE scattering - that the impulse approximation, in which the interactions happen 
only with one nucleon at a time, may not account for the total cross section and that 
instead interactions of the incoming neutrino with nucleon pairs are important. Indeed, 
first results [17, 18, 19, 10] seem to show that the observed high QE scattering cross 
sections could possibly be explained by adding these 2p2h contributions in the primary 
interaction to the impulse approximation cross section. Taking these contributions into 
account would increase the calculated QE-like cross section and thus worsen the agree- 
ment for the 7T/QE ratio. Only when one assumes that a similarly strong contribution 
also to pion production the ratio would stay unaffected. However, such an increase by 
about 25% would still be too low to explain the observed dramatic discrepancy of 60% 
between the GiBUU calculations and the experimental values. 

Another possible explanation for the discrepancy could, of course, be a too low 
neutrino flux assumed in obtaining these pion data, because that would affect both QE 
scattering and pion production. In this view it is interesting to notice, that the agreement 
between our calculations and data is much better for flux-independent and shape-only 
observables, than for the absolute values. 
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